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ABSTRACT 

Context. One of the Friedmann-Lemaitre-Robertson- Walker (FLRW) models that best fits the Wilkinson Microwave 
Anisotropy Probe (WMAP) sky maps of the cosmic microwave background (CMB) is that whose comoving space is 
the Poincare dodecahedral space. The optimal fit of this model to WMAP data was recently found using an optimal 
cross-correlation method. For geometrical reasons, there was concern that systematic error in the estimate of the 
matched-circle (observer-centred) angular radius a, or equivalently, the (comoving) size of the Universe 2rinj (twice the 
injectivity radius), might be much higher than the random error. 

Aims. In order to increase the falsifiability of the model, especially by multiple imaging of collapsed objects, it would 
be useful to reduce the uncertainty in this estimate and to estimate the fraction of the sky where multiply imaged 
gravitationally bound objects should potentially be detectable. 

Methods. A corollary of the matched circles principle — the existence of matched discs — is introduced in order to describe 
a useful subset of multiply imaged objects. The cross-correlation method at ^ 1 h^^ Gpc is applied to WMAP 7-year 
data in order to improve the estimate of a. 

Results. The improved matched-circle radius estimate is a = 23 ±1.4°, where the uncertainty represents systematic error 
dependent on the choices of galactic mask and all-sky map. This is equivalent to 2rinj = 18.2 ±0.5 Gpc for matter 
density parameter ilm ~ 0.28 ± 0.02. The lowest redshift of multiply imaged objects is z = 106 ± 18. Multiply imaged 
high overdensity (rare) peaks visible during 200 > z > 106 should be present in matched discs of radius 14.8 ± 2.3°. 
Conclusions. The accuracy in the matched circle radius estimate is considerably improved by using the higher resolution 
signal. The predicted matched discs (over 200 > z > 106) project to about 20% of the full sky. Since any object located 
exactly in the discs would be multiply imaged at equal redshifts, evolutionary effects would be small for objects that 
are nearly located in the discs. 
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1. Introduction structure on the largest scales is most simply seen in the 

_^ spatial or angular two-point autocorrelation function of the 

I. . The near lack of 10 h Gpc scale structure in cos- temperature fluctuations, as spherical harmonical analysis 

^ ; mic microwave background (CMB) temperature fluctu- requires a "conspiracy" between different harmonics in or- 

^ ation maps was tentatively detected by the COsmic jer to match the obs erved weak correl ation at the largest 

O ■ Background Explorer (COBE) and supported by the scale s (Spe rgel et al. mm Copi et al. \m7[ Sarkar 

' ■ Wilkinson Microwave Anisotropy Probe (WMAP) (Bennett et al. &3TT\ Copi et al. l^DTU)) . 

et al. 12003]). This has long been considered as a sign that The locally homogeneous 3-manifolds that have most 

the Universe, phenomenologically well-approximated as a frequently been studied recently as candidates of comov- 

Friedmann-Lemaitre-Robertson- Walker (FLRW) model, is j^g gpace in order to fit the WMAP data w ith an FLRW 

spatially small. Structures bigger than comoving space it- mod el incl ude the 3-tor us (Sperg el et al. Aurich 

self cannot exist physically, and their observational statis- ^t al. [20071 Aurich 120081 Aurich et al. 120081 l20T0l Aslanyan 

tical description in the universal covering space (ap parent ^ Manohar mWi and the Poincare dode cahedral space 

space) shoul d appr oximately reveal this (StarobinskylMl 53/ j, (Lumi net et al. 1^0051 Auric h et a l. mUM IMbl 

Stevens et al 11993]), even though the apparent nature of the Gundermann 120051 Caillerie et al. 120071 Roukema et al. 

covering space does permit str ucture s that are larger than |2008bl EOPgil) . A full pre- WMAP classification of the pos- 

the Universe (Fig. 1, Roukema UMl) • Sm aU FL RW models jti^ely curved spaces (including S^r) and mathematical 

have long b een discu ssed ( e.g., de Sitt er im ^ Friedmann ^oig foj. studying them in the cosm ological context were 

mm 11924) Lemaitre [19271 Robertson [1935|)|J The lack of presented by Gausmann et al. (1MJT|) . Some analyses sug- 

gest that the Poincare space may not provide a better fit 

^ Lemaitre p927p is the paper in which Lemaitre used redshift than the infinite flat model or a hypersphere modeQ (Key 

and distance estimates of 42 galaxies to estimate the value of — — — — — — — 

what is now known as Ho, finding 575 km/s/Mpc (unweighted) (|2011|l and Block (I2011|l and references therein for discussion of 

or 625 km/s/Mpc (weighted). An incomplete English transla- this issue. 

tion that excludes almost all of the observational analysis was ^ The hypersphere is frequently referred to using the ambigu- 

published four years later (Lemaitre I1931|) . See van den Bergh ous term "closed". 
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I"" matched discs 

•"inj = 7^/10 Rc 

Fig. 1. Relation of matched discs to matched circles for 
the Poincare dodecahedral space, with injectivity radius rinj = 
(7r/10)-Rc, shown in the universal covering space S'' of radius 
_Rc- Multiple copies of the surface of last scattering (SLS) of 
radius rsLS intersect in circles lying in (flat) 2-planes that are 
orthogonal to the plane of the page. The interiors of the circles 
within these 2-planes constitute matched discs. The matched- 
circle (observer-centred) angular radius is a. Redshifts on either 
matched disc increase from Zmin at the centre to zsls on the 
boundary (the matched circle). The SLS can be replaced by 
the 2-sphere defined by z = 200 to obtain matched discs with 

200 > Z > Zmin instead of ZSLS > Z > Zmin- 



et al. 120071 Niarchou & Jaffe I2007p . A heuristic argument 
in favour of the Poincare space is that the residual gravity 
effect implied by imperfect cancellation of the gravitational 
signal from distant images of a nearby massive object sin- 
gles out the Poincare space as being better balanced than 
other spaces (Roukema et al. 120071 Roukema & Rozahski 

mm . 

If the Poincare space hypothesis is correct, then the set 
of density perturbations in the Universe would not consti- 
tute a statistical ensemble of realisations of a model of den- 
sity perturbation statistics, it would consist of a single phys- 
ical realisation within a fundamental domain with a definite 
orientation and comoving size. Thus, it should be possible 
to find the optimal orientation and size of the fundamental 
domain in an astronomical coordinate system, rather than 
only estimate probabilities based on ensemble averages of 
simulated temperature fiuctuations. The parameters of the 
optimal solution should then be usable to design observ- 
ing strategies that should strengthen evidence either for or 
against the hypothesis. 

Estimates of the optimal Poincare space coordinates 
have been made (Roukema et al. I2008bl I2008ap . However, 
these estimates favoured the use of the moderately large- 
scale signal, up to ~ 4 Gpc. For geometrical reasons, 
this had the consequence that one of the fundamental do- 
main parameters, the matched-circle observer-centred an- 
gular radius a, or equivalently, the (comoving) size of the 
Universe (2rinj, i.e. twice the injectivity radius, i.e. the 
length of the shortest closed spatial geodesic), was, in prin- 
ciple, only weakly constrained (Sect. 4.2, Fig. 9, Eq. (32), 
Roukema et al. I2008b[) . The random error in the estimate 
a ~ 21 ± 1° (Table 2, Roukema et al. l2008bp was suspected 
to be dominated by a systematic error of ^ 10°. 

This uncertainty has critical consequences for the falsifi- 
ability of the model. If 32° ^ a <, 48°, then multiply imaged 



objects should be observable at 40^z>,10, respectively, 
near the centres of matched discs (Fig. [IJ. Matched discs 
and this redshift calculation are explained in Sect. 12.11 and 
the results are shown in Fig. [2] This redshift range is ap- 
proximately that when 0.1% or so of baryons are expected 
to have formed H2-cooled protogalaxies within dark mat- 
ter haloes of about 10^ to IO^Mq, in which very high mass 
Population III stars form and some explode as supernovae, 
while others may collapse directly into black holes (e.g.. 
Glover HOnSl Ciardi & Ferrara lMJSl Greif et ah [MM and 
references therein). Although the highest- redshift, spcctro- 
scopically confirmed, collapsed astrophysical object is a 
galaxy at z = 8.6 (Lehnert et al. I2010p . it is generally ex- 
pected that Population III stars may be detected as GRBs 
at 2: > 10 during the coming decades. 

However, if the Universe is larger, i.e. if 18° <, a <, 32°, 
then the corresponding redshifts are higher, in the range 
200 ^ 2: ^ 40, respectively, so protogalaxies can have col- 
lapsed only in higher overdensity peaks, i.e. much more 
rarely. Their rareness at higher redshifts decreases the 
chance that observable objects will lie near the centres of 
the matched discs compared to the case for larger values of 
a. 

The advantage of seeing two topologically lensed images 
of an object at nearly the same epoch is that evolutionary 
effects would be minimal, although there would still be a 
projection effect induced by seeing the same object "from 
front and back". Significant evidence of the absence of a 
predicted topologically lensed image would be able to re- 
fute the Poincare dodecahedral space model. Alternatively, 
well-matched topologically lensed pairs would increase the 
over-determined nature of the model, leading to succes- 
sively more and more topological lensing predictions. Thus, 
it is important to see if the estimate of a can be made more 
accurate. 

A more accurate estimate of a is made by using the 
cross-correlation method (Roukema et al. I2008bl I2008al) at 
higher comoving spatial resolution ( <, 1 h^^ Gpc) than in 
the earlier work and updating to the WMAP 7- year data re- 
lease. The method is presented in Sects l2.2l and l2.3l Results 
are presented in Sect.[2]and conclusions are given in Sect. 21 
All distances are FLRW comoving distances. The Hubble 
constant is written Hq — lOO/i km/s/Mpc. 

2. Method 

2.1. Matched discs 

Here we introduce an observational corollary of the iden- 
tified circles principle (Cornish et al. 119961 119981 see also 
Fig. 2, Roukema I2000p . Two circles on the surface of last 
scattering (SLS) are identified because in the comoving cov- 
ering 3-space, they consist of the intersection of two copies 
of a single physical (fiat) 2-plane with a 2-sphere (the SLS). 
Filled polygonal subsets of the two copies of the 2-plane 
constitute two faces of a copy of the fundamental domain 
that can be thought of as being pasted together. The flat 
interiors of the two matching circles — within the two copies 
of the flat 2-plane — constitute matched discs with sub-SLS 
redshifts, i.e. Zmin l£ z < zg^s, where Zmin occurs at the 
disc centre. This is illustrated in Fig. [TJ 

A large enough a (small enough 2ri„j) would im- 
ply a low enough z^in for collapsed astrophysical objects 
[e.g. quasars, Lyman-break objects or gamma-ray bursts 
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Fig. 2. Upper panel: Redshift Zmin at the centre of a matched 
disc in a Poincare space as a function of matched circle radius 
a, for = 0.26, 0.28, 0.30, 0.32 as curved sohd, dotted, dot- 
dashed and dashed hues, respectively (almost indistinguishable). 
Redshifts Zmi-n ~ 10, 40, and 200 are indicated as illustrations of 
epochs of successively rarer expected formation of galaxies with 
Population III stars. Lower panel: Corresponding injectivity ra- 
dius of the Universe rinj . 



(GRBs) from early forming galaxies] located near the cen- 
tres of these matched discs to be seen as nearly antipodal 
pairs at nearly identical redshifts. Objects successively fur- 
ther from the centres of the matched discs would be multi- 
ply imaged at successively higher redshifts z > Zmin, i-e. at 
successively earlier, but still approximately equal, epochs. A 
survey whose selection criteria select a redshift band beyond 
Zmin should detect an annulus of multiply imaged objects. 
For example, if z^in = 20 and the survey selects objects in 
the range 30 < 2; < 40, then there should be an annulus of 
multiply imaged objects. 

The redshift to the centres of a pair of matched discs for 
a given matched circle radius a can be calculated by find- 
ing the value of the dark energy parameter Oa that equates 
the curvature radius calculated locally with its value cal- 
culated globally using the injectivity radius. The former is 
Rc = (c/Hq) (17ni + JIa — where c is the spacetime 

unit conversion constant, is a fixed value of the matter 



density parameter, and Hq is the Hubble constant. The cal- 
culation using the injectivity radius is made using Eq. (15) 
of Roukema et al. (|2008bp . a and the comoving distance 
to the surface of last scattering rsLs(f^m, ^A, Hq, zsls), as- 
suming zsLS = 1100. Thus, 



Rc = 



1 



Hq + ri,„ — 1 





11,^ A, Hq) 


arctan 


"tan(7r/10)" 
cos a 



(1) 



Given fim and Hq, the latter equality can be solved numer- 
ically for yielding i?c from the former equality. The 
injectivity radius is 



= (^/10)il'c 



(2) 



so the redshift 2,„in to the centre of the matched discs is 
obtained by inverting 



r(r2ni,f^A,-ffo,2:min) 



(3) 



where r is the standard FLRW radial comoving distanceO 
Figure[2] (upper panel) shows that if 32° ^ a ^ 48°, then 
the redshifts of objects lying in the matched discs would be 
in the range 40 z >, 10, respectively. For smaller matched 
circle radii 18°;$a;S32°, the corresponding redshifts are 



higher, in the range 200 , 



,40, respectively. 



Larger matched circle radii correspond to lower red- 
shifts, but also higher total density and smaller curvature 
radii Rc- For example. Fig. 15 of Gausmann et al. (|2001l) 
shows a multiple imaging spike with a minimum redshift 
lying in the range 1 < z < 3 when fim = 0.35, Ha = 0.75. 
This can be calculated by using the left-hand equality 
of Eq. (P to obtain Rq = 9A8h-^ Gpc, Eq. (0) to ob- 
tain Tinj = 2.98/i~^ Gpc, and inverting Eq. ([3]), yielding 
z = 1.45. Use of the right-hand expression in Eq. ([T|) shows 
that the corresponding matched circle radius is a = 77°. 

The lower panel of Fig.[2]shows the corresponding injec- 
tivity radius rinj of the Universe. While the redshift at the 
centre of a matched disc is insensitive to fim , the size of the 
Universe is sensitive to fim- When a = 0°, the injectivity 
radius is identical to the comoving distance to the surface 
of last scattering, i.e. ri„j = rsLS [cf. Eqs ([1]) and 



2.2. Cross-correlation maximisation method 

To improve the estimate of a, the cross-correlation method 
(Roukema et al. I2008b[ I2008ap is used over the single pa- 
rameter a, using the previously found solution for the po- 
sitions of the axes of the fundamental domain. The cross- 
correlation method can be thought of as reconstructing the 
spatial two-point correlation function of density fluctua- 
tions (using temperature fluctuations as a proxy) at sub- 
injectivity radius length scales, where an observed pair of 
points is only used if the observed spatial separation in 
the comoving covering space is larger than the injectivity 
radius. Holonomies gi are applied to one member ?/ of a 
pair of observed spatial points {x,y) so that if the implied 
pair {x,g~ (y)) is a close pair, then the observed temper- 
ature fluctuation product AT{x) AT{y) is cumulated to 
the cross-correlation estimate at the separation d[x, g~^ (y)] . 
The method can also be thought of as an extension of the 
identified circles principle (Cornish et al. 119961 11998p . but 
with a uniform sky weighting over observed points instead 



It might be useful to call Zmin the injectivity redshift Zinj. 
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Fig. 3. Gigaparsec-scale cross-correlation [Eq. ((6|] in 

^K^ as a function of matched-circle radius a, for the WMAP 
7-year ILC map, for fractions of the surface of last scattering 
radius /3 — 0.03, 0.1, 0.2, 0.4 as indicated (from top to bottom at 
a ~ 23°), using the KQ85 galactic mask. 
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Fig. 4. Gigaparsec-scale cross-correlation £,0{a), as for Fig. [31 
for the WMAP7 foreground-corrected V map. 




Fig. 5. Gigaparsec-scale cross-correlation ^/3(a), as for Fig. [3l 
for the WMAP7 foreground-corrected W map. 
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Fig. 6. Gigaparsec-scale cross-correlation ^/3(a), as for Fig. [3l 
for the WMAP7 ILC map and the KQ75 galactic mask. 



of a weighting that favours matched circle intersections. 
If a hypothesised set of sky orientation, injectivity radius 
and twist parameters are physically correct, then the lo- 
cally isotropic components of the auto-correlation function 
should be reproduced by the cross-correlation function. 

2.3. Correlation functions and optimisation criterion 

The Markov chain Monte Carlo ( MCMC) method used 
in Roukema et al. (|2008bl I2008al) adopted Np = 2000 
points per position in MCMC parameter space (Sect. 3.5, 
Roukema et al. I2008bp . The optimisation criterion (a 
pseudo-probability) was aimed at maximising — 
at <,4 Gpc scales. The definition of the criterion 
[Eqs (25), (26), Roukema et al. I2008b| took into account 
the uncertainty in estimating the cross-correlation (f) in- 
duced by using a small value of Np . This provided a practi- 
cal method of exploring a large volume of parameter space, 
focussing on cross- and auto-correlations mostly at the up- 
per end of the separation range r <, 4 h^^ Gpc. 

Here, in order to optimise the estimate of a, a much 
higher number of points is used per correlation calculation, 
i.e. A^p = 2 X 10^, and only a is varied. This enables the 
use of the correlation signal at separations r <, 1 Gpc. 



From Fig. 11 of Roukema et al. (j2008b[) . it is clear that 
^ 'fc at r <, 1 h^^ Gpc. Because of the complicated mix 
of different effects — the naive Sachs- Wolfe effect, the inte- 
grated Sachs- Wolfe effect, the Doppler effect, and intrinsic 
temperature fluctuations — that contribute to and the 
auto-correlation ^a, it is not obvious whether should be 
as high as ^a at sub-Gpc scales. In other words, on these 
smaller scales, the increasing contribution of the Doppler 
effect — which is locally anisotropic — to the temperature 
fluctuations makes it less realistic to test the hypothesis 
that and ^a are identical functions apart from calcula- 
tional uncertainty. Moreover, the higher value of A^p reduces 
the numerical usefulness of estimating the calculational un- 
certainty. 

Hence, with the aim of maximising the cross-correlation 
of locally isotropic emission, the optimisation criterion 
is defined as the mean cross-correlation below a given 
length scale, without a direct comparison with the auto- 
correlation. However, since the comoving (spatial geodesic) 
separation d{9A, ^'sls) between two fixed points on the SLS 
separated by an observer-centred angle 9a varies with Jim 
and SI A, d needs to be replaced by a scaled separation /3 
that to a good approximation depends only on 9a, in or- 
der that Ca(/3) has only a weak dependence on r2,„ and f^A- 
This can be done for small separations d, which are the ones 
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800 




Fig. 7. Spatial auto- and cross-correlation functions, and 
5c, respectively, for the WMAP7 ILC map for matched-circle 
radius a = 23°, shown as thin and thick (colour online) curves, 
respectively. The cross-correlation function is uniformly sampled 
only for /3 < 0.44 (for calculation speed); at higher separations, 
it is, in principle, biased. In order to better compare the small- 
scale behaviour of and , the auto-correlation reduced (arbi- 
trarily) by a factor of 4 is also shown. The mostly upper (lower) 
member of each close pair of functions is for the KQ85 (KQ75) 
galactic mask. 



of most interest. As in Eq. (16) of Roukema et al. (j2008bp . 
the spherical sine law gives 



sin[d(0A,rsLs)/(2i?c)] 



rshS 



1 . 0A . J'SLS 

sm — sm , 

'^SLS 2 Rc 



I.e. 



d(6'A,r-sLs) 



rsLS 



„ -Rc . »'SLS . 0A 

2 sm— — sm— , 

^SLS -Kc ^ 



(4) 



(5) 



the 



where the second and higher order terms in 
Taylor expansion of sin[(i/(2i?c)] are dropped, since 
for d<lh-^ Gpc and Rg > 25 h- \ Gpc, d/{2Rc) < 0.02. 
Using Eq. (15) of Roukema et al. (j2008bl) . the ratio ^ 
depends only on a, without any direct dependence on ilm 
and Qa- For 0° < a < 45°, we have 1.967 > 2-^ sin ^ > 
1.939, respectively, i.e. the latter expression varies by less 
than 1.5% over this range. In contrast, for 0.25 < flm < 
0.35, rsLS varies by about 14%. Thus, /3 := d/rsLS depends 
mostly on 9a and only weakly on fl^ and flA- 

This is used for the mean cross-correlation, i.e. the lat- 
ter is calculated below a fixed fraction /3 of the comoving 
distance to the surface of last scattering radius, 



1 



(6) 



The primary scales of interest are (3 = 0.033 and 0.1, cor- 
responding to /3rsLS ~ 0.33 and 1 Gpc, respectively. 

The primary estimate is the matched circle radius a 
that maximises ^i3{a) for /? = 0.033 and 0.1 for the 7- year 
WMAP ILC map|3 masking with the KQ85 mask to min- 
imise galactic contamination. Estimates of the systematic 
error, assumed to dominate the random error, are made by 



|http : //la mbda . gs f c . nasa . gov/dat a/map/dr4/ dip/ \ 
|ilc/vnnap_ilc_7yr_v4. f its| 



Table 1. Matched-circle radius a at which the Gpc-scale cross- 
correlation 5/3 (q) is maximised. 



map 




/3 








0.033 


0.1 


0.2 


0.4 


ILC, 85" 


23.0 


23.0 


23.0 


23.0 


V, 85" 


21.8 


23.0 


24.4 


24.4 


W, 85" 


22.0 


22.4 


22.4 


34.4 


ILC, 75** 


23.8 


23.8 


23.8 


23.8 



Notes. ("' KQ85 mask KQ75 mask 

considering the foreground-reduced V banc0 and W band 
mapslfl and by applying the KQ75 mask. 

Since the aim here is to see whether the uncertainty 
in a can be reduced so that it no longer dominates the 
uncertainty in the optimal solution, the sky orientation pa- 
rameters of the fundamental domain are fixed at (l, b, 6) — 
(184.0°, 62.0°, 34.0°), usin g the ex act dodecahedral solution 
found in Roukema et al. (j2008bp (Fig. 11 caption) rather 
than the individual axis estimates. 

3. Results 

3.1. Matched circle radius a 

Figure P shows a well-defined local peak in the mean corre- 
lation ^p{a) near a = 23°. Even for /3 = 0.4, i.e. using the 
<Ah^^ Gpc cross-correlations, a small local peak near the 
same value is present. The smaller scale cross-correlations 
clearly overcome the concern of a systematic error of about 
10° in the estimate of a. 

Several features of this figure can be interpreted sim- 
ply if the overall Poincare dodecahedral space solution is 
physically correct. Firstly, as would be expected for a true 
solution, the rough solution found using the larger scale sig- 
nal becomes better defined using the smaller scale signal. 
Secondly, as a decreases down to 10°, f/3(a) drops strongly. 
The smaller [3 is, the sharper the drop. This is consistent 
with the solution being correct, since matched circles for 
wrong values of a that are smaller than the would-be true 
value force small-scale cross-correlations to be estimated at 
large true separations, where the auto-correlation is weak. 

Figures |4] and [5] show that the peak in (a) is present 
in both the V and W foreground-reduced maps. Figure [5] 
shows that the peak is a little weakened using a stronger 
galactic cut, but remains present. Because of the 36° rota- 
tion used in matching circles or discs, stronger galactic cuts 
severely reduce the numbers of matched pairs and nearly 
matched pairs. 

The matched circle radii that maximise ^p, after a 
Gaussian smoothing with ctq ~ 0.2°, are listed in Table[TJ In 
the ILC KQ85-masked map, the maximum cross-correlation 
is independent of the integration length scale /3, giving a — 
23.0° in all four cases. A sharp maximum cross-correlation 
is clearly present at about the same matched circle radius 
for the single-band V and W (foreground-reduced) maps, 
but with some dependence on /3. The ^ 1 Gpc scale 



^ http : / /lambda. gsf c .nasa. gov/dat a/map/dr4/ skymaps/ 
7yr / f orer ed/ wmap_band_f orered_iqumap_r9_7yr_V_v4 .fits 

" |h ttp : / /lambda. gsf c .nasa. gov/dat a/map/dr4/ skymaps/J 
|7yr/f orered/wmap_band_f orered_iqumap_r9-7yr_W_v4. f its, 
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V and W maps (/3 = 0.033,0.1) give 21.8° < a <_23°. 
The ILC KQ75-masked map again gives a maximum £,i3{a) 
that is independent of /3, i.e. a = 23.8°, a little higher than 
the other estimates. Conservatively taking the maximum 
differences to be independent and adding them in quadra- 
ture, the systematic error in the estimate of a induced by 
the choice of the map and the galactic mask is 1.4°, giving 
a = 23.0 ± 1.4° as the improved estimate of the matched 
circle radius. 

The cross-correlation functions for the ILC map for 
a = 23° are shown in Fig. [71 This should be compared 
with Fig. 11 of Roukema et al. (|2008bp . In Fig. El since 
the method of Roukema et al. (|2008ap is used to speed up 
the calculations, only separations below (3 = 0.44, i.e. about 
4.4 Gpc, are based on an underlying uniform sky distri- 
bution. Nevertheless, cross-correlations at higher scales are 
similar to those in Fig. 11 of Roukema et al. (I2008bp . i.e. 
they do not appear to be strongly affected by the potential 
bias. 

The main new feature in Fig. [71 is that there is a sharp 
upturn at sub-Gpc scales, which was not visible in Fig. 11 
of Roukema et al. (|2008bp . rather than a continuation of 
the approximately constant slope on 1-5 Gpc scales, 
visible in both figures. If the naive Sachs- Wolfe effect, i.e. 
a locally isotropic effect, contributes about half of typical 
temperature fluctuation amplitudes on these scales, then 
the expected cross-correlation should be ^ Ca/4. The 
latter function is plotted in Fig. [71 in addition to ^a- This 
shows that Ca/4 indeed gives a rough estimate of the 
cross-correlation on sub-Gpc scales with /3 > 0.011, i.e., 
firshS ~ 0.11 Gpc. Thus, a possible order-of-magnitude 
interpretation would be that about half the temperature 
fluctuations on these scales are locally isotropic. 

3.2. Inferred parameters: 2rinj, Zmin, q;2oo 

As described above, for a choice of fim, H^, and a, Eqs ^ 
and ([2]) can be used to flnd SIa, i?c, Hnj and the redshift 
to the matched disc centres Zmin- For Q,^ = 0.28 ±0.02 and 
a = 23 ± 1.4°, and assuming that these give Gaussian inde- 
pendent errors, twice the injectivity radius of the Universe 
is 2rinj — 18.2 ± 0.5 h^^ Gpc, and the redshift to the 
matched disc centres is Zmin = 106 ± 18. 

The angular radius of a matched disc from its centre 
at Zrain = 106 ± 18 up to a redshift at which some of the 
very high density peaks should first collapse, z = 200, can 
be calculated analogously to Eq. (15) of Roukema et al. 
mom , i.e. 

tan(rsLs/-Rc)l 

q;9oo = arccos cos a -. tttt i (') 

[ tan(r2oo/i?c) J 

where r2oo is the comoving radial distance to z = 200 for the 
same values of rim, ^^a, and Hq as for the other parameters. 
This yields a2oo — 14.8 ± 2.3°. These discs do not intersect 
with one another. Thus, for the full sky, the 12 matched 
discs for 200 > z > 106 would project to a fraction of 
247r(l-cosa20o)/(47r) = 20 ±6% of the sky. The fractional 
coverage of the unmasked sky should be similar. 

4. Conclusion 

The sub-Gpc cross-correlation estimates shown in Figs [SHU 
show that the systematic uncertainty in the earlier esti- 
mate of the matched circle radius for the optimal fit of the 



Poincare dodecahedral space model to WMAP sky maps 
is signiflcantly reduced by using this high-resolution signal, 
giving a = 23 ± 1.4°. For = 0.28 ± 0.02, the inferred 
size of the Universe is 2rinj = 18.2 ± 0.5 Gpc. 

This constraint on a pushes sub-SLS multiple-image 
testing of the Poincare space model (at least) several 
decades into the future, given present plans for new tele- 
scope projects. Nevertheless, higher density peaks should 
collapse and form population III stars earlier than lower 
density peaks, so future surveys of very high density peaks 
that collapse early, during 200 > z > 106, should, in prin- 
ciple, be feasible. The objects would be (multiply) visible 
in 12 matched flat discs that project to about 20% of the 
full sky. The redshifts of the objects increase radially from 
z = 106 ± 18 at the centre of a disc out to the edge of 
the disc defined by z = 200. An annulus extending the 
disc further to zsls ~ 1100 on the SLS would define a 
multiply imaged region where collapsed objects would be 
extremely rare. The two topologically lensed images of a 
physical object are located at identical redshifts at widely 
separated directions in the sky in discs that are antipodal 
and matched with a twist of -|-36° (right-handed). The red- 
shift equality of the two images of an object implies that 
evolutionary effects would not be a problem for objects lo- 
cated exactly in the discs. Objects that are in the close fore- 
ground/background (foreground in one disc, background in 
the other disc) would have small evolutionary differences 
between the two images of an object. 

Acknowledgements. A part of this project has made use of Program 
Obhczen Wielkich Wyzwah nauki i techniki (POWIEW) compu- 
tational resources (grant 87) at the Poznah Supercomputing and 
Networking Center (PCSS). Use was made of the Centre de Donnees 
astronomiques de Strasbourg (http://cdsads.u-strasbg.fr), 
the GNU PLOTUTILS graphics package, and the GNU OCTAVE 
command-line, high-level numerical computation software 
( http : / /www . gnu . org/sof tware/octave ) . 



References 

Aslanyan, G. , &c M anohar, A. V. 2011, ArXiv e-prints, 

[arXiv: 110 4.0015] I 
Aurich, R. 2008, ClassQuantGra, 25, 225017, [arXiv: 0803. 2130] 
Aurich, R., Janzer, H. S., Lustig, S., & Steiner, F. 2008, Classical and 

Quantum Gravity, 25, 125006, [arXiv : 0708 . 1420] 
Aurich, R., Lustig, S., & Steiner, F. 2005a, ClassQuantGra, 22, 3443, 

[arXiv : astro-ph/0504656] 
Aurich, R., Lustig, S., & Steiner, F. 2005b, ClassQuantGra, 22, 2061, 

[arXiv : astro-ph/ 0412569] 
Aurich, R., Lustig, S., & Steiner, F. 2010, ClassQuantGra, 27, 095009, 

[arXiv: 0903. 3133] 
Aurich, R., Lustig, S., Stei ner, F., fc The n, H. 2007, ClassQuantGra, 

24, 1879, [arXlv:astro-ph/0612308] 
Bennett, C. L., Halpern, M., Hinshaw, G., et al. 2003, ApJS, 148, 1, 

[arXiv:astro-ph/0302207] 

Block, D. L. 2011, ArXiv e-prints, [ arXiv: 1106. 3928] | 

Caillerie, S., Lachicze-Rcy , M., Luminet, J.-P., et al. 2007, A&A, 476, 

691, [arXiv:0705.021 7v2]1 
Ciardi, B., & Ferrara, A. 2005, Space Sci. Rev., 116, 625, 

[arXiv : astro-ph/04090 18] 
Copi, G. J., Huterer, D., Schwarz, D. J., & Starkman, G. D. 2007, 

Phys. Rev. D, 75, 023507, [arXiv : astro-ph/0605135] 
Copi, C. J., Huterer, D., Schwarz, D. J., & Starkman, G. D. 2009, 

MNRAS, 399, 295, [arXiv: 0808. 3767] 
Copi, C. J., Huterer, D., Schwarz, D. J., & Starkman, G. D. 2010, 

Advances in Astronomy, 2010, 847541, [arXiv: 1004.5602] 
Cornish, N. J., Spergel, D. N., & Starkman, G. D. 1996, ArXiv 

Gen.Rel. & Quant. Gosm. e-prints, [arXiv : gr-qc/9602039] 
Cornish, N. J., Spergel, D. N., & Starkman, G. D. 1998, 

ClassQuantGra, 15, 2657 
de Sitter, W. 1917, MNRAS, 78, 3 



6 



Roukema & Kazimierczak: Poincare space Universe size 



Priedmann, A. 1923, Mir kak prostranstvo i vremya (The Universe as 

Space and Time) (Leningrad: Academia) 
Friedmann, A. 1924, Zeitschr. fiir Phys., 21, 326 

Gausmann, E., Lelioucq, R., Luminet, J. -P., Uzan, J. -P., & Weeks, J. 

2001, ClassQuantGra, 18, 5155, [arXiv : gr-qc/0106033] 
Glover, S. 2005, Space Sci. Rev., 117, 445, [arXiv:astro-ph/0409737] " 
Greif, T. H., Johnson, J. L., Klessen, R. S., & Bromm, V. 2008, 

MNRAS, 387, 1021, [arXiv: 0803. 2237] 
Gundermann, J. 2005, ArXiv e-prints, [arXivT astro-ph/0 503014] | 
Key, J. S., Cornish, N. J., Spergel, D. N., & Starkman, G. D. 2007, 

Phys. Rev. D, 75, 084034, [arXiv:astro-ph/0604616] 
Lehnert, M. D., Nesvadba, N. P. H., Cuby, J.-G., ct al. 2010, Nature, 

467, 940, [arXiv: 1010.4312] 
Lemaitre, G. 1927, Annales de la Societe Scientifique de Bruxelles, 47, 

49 

Lemaitre, G. 1931, MNRAS, 91, 483 

Luminet, J., Weeks, J. R., Riazuelo, A., Lehoucq, R., &; Uzan, J. 2003, 

Nature, 425, 593, [arXiv: astro-ph/0310253] , 
Niarchou, A., fc Jaffe, ^OOT, FRL, 99, 081302, 

[arXlv : as tro-ph/0702436]'' 
Robertson, H. P. 1935, ApJ, 82, 284 

Roukema , B. F. 1996, MNRAS, 283, 1147, 

[arXiv : astro-ph/9603052] 
Roukema, F. 2000, Bull. Astr. Soc. India, 28, 483, 

[arXiv : astro-ph/0010185] 
Roukema, B. F., Bajtlik, S., Biesiada, M., Szanicwska, A., & 

Jurkiewicz, H. 2007, A&A, 463, 861, [arXiv:astro-ph/0602159] 
Roukema, B. F., Buliriski, Z., & Gaudin, N. E. 2008a, A&A, 492, 673, 

[arXiv: 0807. 4260] 
Roukema, B. F., Bulinski, Z., Szanicwska, A., & Gaudin, N. E. 2008b, 

A&A, 486, 55, [arXiv: 0801. 0006] 
Ro ukema, B. F., & Roiahski, P. T. 2009, A&A, 502, 27, 

I [arXlv"T09 02 . 3402] 
Sarkar, D., Huterer, D., Copi, C. J., Starkman, G. D., & Schwarz, 

D. J. 2011, Astroparticle Physics, 34, 591, [arXiv: 1004.3784] 
Spergel, D. N., Verde, L., Pciris, H. V., et al. 2003, ApJS, 148, 175, 



I [arXiv:astro-ph/0302209] 
Starobinsky, A. A. 1993, JETPLett, 57, 622 

Stevens, D., Scott, D., & Silk, J. 1993, PRL, 71, 20 

van den Bergh, S. 2011, JRASG, in press. . [arXiv : 1106 .1195] I 



7 



